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ABSTRACT 

We introduce a method for calculating the probability density function (PDF) of a 
turbulent density field in three dimensions using only information contained in the 
projected two-dimensional column density field. We test the method by applying it 
to numerical simulations of hydrodynamic and magnetohydrodynamic turbulence in 
molecular clouds. To a good approximation, the PDF of log (normalised column den- 
sity) is a compressed, shifted version of the PDF of log (normalised density). The 
degree of compression can be determined observationally from the column density 
power spectrum, under the assumption of statistical isotropy of the turbulence. 
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1 INTRODUCTION 

The probability density function (PDF) of the density field 
in molecular clouds is a key ingredient in most analytic mod- 
els of star formation (Padoan & Nordlund 2002; Krumholz & 
McKee 2005; Elmegreen 2008; Hennebelle & Chabrier 2008; 
Padoan & Nordlund 2009; Hennebelle & Chabrier 2009). 
Knowledge of the density PDF is required to determine the 
overall star formation rate or efficiency. In some models, the 
density PDF is of central importance in determining the 
emergent stellar initial mass function. The majority of mod- 
els assume a lognormal density PDF (Vazquez-Semadeni 
1994), with the width of the PDF increasing with the Mach 
number of the turbulence (Padoan, Nordlund, & Jones 1997; 
Passot & Vazquez-Semadeni 1998; Federrath, Klessen, & 
Schmidt 2008). 

Observational knowledge of density fields in molecular 
clouds (let alone the PDF) is very limited. We do not have 
access to the density field in three dimensions (3D) but in- 
stead can only view the projected column density field in two 
dimensions (2D). Use of molecular tracers of different critical 
density can in principle yield some information, but even the 
most sophisticated excitation analyses provide only a single 
"density" per line-of-sight whereas the transverse variations 
in "density" so measured must imply comparable, or per- 
haps greater, fluctuations in density along the line-of-sight. 
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A route to the PDF could involve (for example) the mea- 
surement of mass exceeding a range of critical densities from 
a suite of tracers. While the amount of data involved here 
would likely be prohibitive for nearby clouds of large angular 
extent, it can be usefully applied in an extragalactic context 
(e.g. Krumholz & Thompson 2007). 

Column density fields traced by extinction of back- 
ground stars are perhaps the most robust way of acquiring 
constraining data on the density PDF in nearby clouds. Col- 
umn density PDFs can be constructed (e.g. Cambresy 1999; 
Lombardi 2009; Kainulainen et al 2009) but the relation be- 
tween these and the 3D density PDF is currently unknown. 
Compression of the PDF due to line-of-sight averaging is ex- 
pected, and there are some indications that a lognormal den- 
sity PDF will project into a (less broad) lognormal column 
density PDF (Ostriker, Stone, & Gammie 2001; Vazquez- 
Semadeni & Garcia 2001, Federrath et al 2009). The degree 
of compression is presently unknown, but recent work by 
Brunt, Federrath, & Price (2010; hereafter BFP) demon- 
strated how to calculate the normalised density variance in 
3D from information contained solely in the column density 
field. Comparison of the measured normalised column den- 
sity variance with the inferred normalised density variance 
can provide some information on the degree of compression. 

In this paper, we introduce a method by which the 3D 
density PDF can be constructed from measurements made 
on the column density field alone. Measurements of the nor- 
malised column density variance, the column density power 
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spectrum, and the column density PDF are required, and 
can be combined to construct an estimate of the 3D density 
PDF. 



2 THE METHOD 

2.1 Reconstructing the 3D Density Variance 

We define the normalised density field in 3D, X3, as: 

X3 = P/PO, (1) 

where p is the density and po is the mean density. Similarly, 
in 2D, the normalised column density, X2, is defined as: 



x 2 = N/No, 



(2) 



where N is the column density and No is the mean column 
density. 

The normalised density variance, a^ 3 is given by: 

°l 3 = (2:3) - {x-i} 2 , (3) 

and the normalised column density variance, a^ 2 , is given 
by: 



(x 2 ) - (2:2} 



(4) 



where angle brackets denote averaging over the fields. 

While the field X3 is observationally inaccessible, BFP 
showed that a^ 3 can be estimated solely from measurements 
made on the field X2- Making use of Parseval's Theorem, 
BFP define the 2D-to-3D variance ratio, R, as: 

/ A/2 A/2 \ 

E E {p*.){k)\ - p*m 

V k^--X/2+l fc H = -A/2+l / 



R = 



A/2 A/2 A/2 \ 

IEEE ( p ^)(k)\ ~p X2 (o) 

yfc x = -A/2+l fe M = -A/2+l k z = -\/2+l I 

(5) 

where {P X2 )(k) is the azimuthally averaged power spectrum 
of X2, P X2 (0) is the power spectrum of X2 evaluated at k = 
0, and A is the scale ratio of the field (i.e. the number of 
pixels along each axis, with the requirement that the field is 
square). 

The inferred 3D density variance, o"^ 3R , can then be 
calculated via: 



and BFP show that: 



(6) 



oizR^ois (7) 

to about 10% accuracy as long as the field X3 is statisti- 
cally isotropic (see BFP for a discussion of this require- 
ment). Straightforward modifications of equation (0) and 
equation © to account for the effect of a telescope beam 
and to account for the effect of zero-padding (to produce 
a square field and/or to reduce edge discontinuities in the 
power spectrum calculation) are given in BFP. 

It is important to recognise that the variance of a field 
calculated at finite resolution is necessarily a lower limit to 
the variance that would be observed in the continuous limit 
(i.e. arbitrarily high resolution). Brunt (2009) discusses the 
consequences of this when the BFP method is applied to the 



Taurus Molecular Cloud, finding that <J X3R may underesti- 
mate the variance that would be obtained in the continuous 
limit by as much as a factor of 2. For the numerical simula- 
tions used in this paper, there is no sub-resolution structure, 
and the 3D variances are not subject to this problem. 

2.2 Reconstructing the 3D Density PDF 

In the following analysis, we propose to make a transforma- 
tion of the field X2 such that the normalised variance of the 
transformed field is equal to <? 2 3R . Our conjecture is that the 
transformed field has the same PDF as the 3D normalised 
density field X3. Note that a suitable transformation can 
match the first two moments of the normalised transformed 
field to the first two moments of the normalised 3D density 
field (mean = unity, variance = (J^ 3R ~ &x 3 )- We do not 
have access to the higher order 3D moments, so must rely 
on experiment to investigate the reliability of this procedure. 

We begin by noting that a transformation X2 — > ax2, 
where a is a constant, has no effect on the normalised vari- 
ance due to the re-normalisation of ax2 ■ The simplest useful 
transformation is therefore X2 — > ax\ where a and £ are 
constants. We therefore define the transformed field, x^r, 
via: 



X3R = ax 2, 



(8) 



where the re-normalising constant, a, is given by: 

a=(x 2 )-\ (9) 

and we have called the transformed field X3R (even though 
it is technically a 2D field) since we intend to endow it with 
a variance of cr 2 3R ■ This can be simply achieved by refining 

a series of test values of £ until the variance of ax\ matches 

2 

If both the density PDF and column density PDF arc 
lognormal in form, the appropriate value of £ can be simply 
derived. First, we note that equation (JSJ) is equivalent to: 



ln(a;3ij) = ln(a) +£ln(x 2 ), 



(10) 



and therefore that the variances of ln(x3n) and ln(a;2) are 
related by: 



(11) 



independently of the value of a. 

With lognormal PDFs for density and column density: 



°~in(x 2 ) = hi(l + o X2 ), 
and 

a ln(x 3R ) = ln(l + °l 3R ) 



so that: 



Mi + 4, J 

ln(l + a| 2 ) 



(12) 



(13) 



(14) 



In general, the PDFs will not be lognormal in form and the 
more general procedure for deriving £ should be employed. 

An important component of the above reasoning (most 
clearly seen in equation (|10|l which applies independently of 
the forms of the PDFs) is that we are assuming that Ixi(xsr) 
is a scaled, shifted version of ln(x2), and therefore that the 
PDF of ln(a;3ij) is a scaled, shifted version of the PDF of 
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\n(x2)- A test of whether the 3D normalised density PDF is 
recoverable by the above method is therefore also a test of 
whether the form of the PDF of ln(a;3) is preserved during 
the projection to 2D. Clearly this will not be true in general 
but may be true of a restricted set of PDFs that characterise 
interstellar density fields. Some support for this is already 
in the literature (Ostriker et al 2001; Vazquez- Semadeni & 
Garcia 2001, Federrath et al 2009). In particular, Vazquez- 
Semadeni & Garcia (2001) suggested that column density 
PDFs could appear lognormal, but tended to appear nor- 
mal if the correlation length of the density field is small 
compared to the line-of-sight extent of the cloud. This is con- 
sistent with our picture: the correlation length is encoded in 
the power spectrum (e.g. by a turnover at some spatial fre- 
quency). The line-of-sight compression of the variance is con- 
trolled by the power spectrum and is more pronounced for a 
flatter spectral slope. If the variance is reduced significantly, 
the density PDF can project into an (apparently) normal 
column density PDF, since a lognormal with very small vari- 
ance is approximately normal. To constrain changes in the 
form of the PDF under projection would require knowledge 
of how additional higher order moments (beyond the second 
given by BFP) are related in 3D and 2D. 

Finally, we note that scaling the PDF of ln(a;3) to the 
PDF of ln(p) = ln(x3po) = 111(2:3) + ln(po) is trivial, as long 
as an estimate of po is available. Lack of knowledge of po 
does not affect the above analysis. However, it is important 
to recognise that a PDF obtained at finite resolution may 
be a compressed version of the true PDF that would be 
obtained in the continuous limit. 



3 APPLICATION TO NUMERICAL 
SIMULATIONS OF TURBULENCE 

3.1 Magnetohydrodynamic Simulations 

We now test the method described above on numerically- 
simulated magnetohydrodynamic (MHD) turbulent density 
fields, obtained at a range of Mach numbers (M = 2 to 20) 
and magnetic field strengths (ratio of thermal-to-magnetic 
pressure f3 =0.00125 to 10). The MHD simulations were 
computed using the grid-based code FLASH (Fryxell et al 
2000) using solenoidal forcing of turbulence at large scales 
(Federrath et al 2009; Price & Federrath 2010; Brunt 2003; 
Brunt, Heyer, & Mac Low 2009). Previously, BFP have 
analysed these fields with the goal of determining the 3D 
normalised density variance using only 2D column density 
maps, so that values of R (equation are available for all 
fields. Here, we use a representative subset of the simulated 
density fields to test the 3D PDF reconstruction method. 

The basic procedure is as follows. Each 3D density field 
(scale ratio A = 256) is integrated over one of the cardinal di- 
rections to produce a 2D column density field, which is sub- 
sequently normalised to produce a field x%, of mean unity. 
The normalised column density variance, 0% , is calculated 
for each field. The power spectrum of X2 is measured, and 
used to calculate R via equation ((5|, which then allows cal- 
culation of o~x 3R = o~i 2 /R. Since the PDFs of these fields are 
approximately lognormal, we calculated £ via equation (| 14p 
and then transformed X2 into X3H via equation (SJ. Finally, 
the PDF of ln(a;3_R) is calculated for comparison to the PDF 



of ln(a;3) where 2:3 is the true normalised 3D density field. 
The PDFs of ln(x3fi) and ln(a;3) are calculated over the same 
range and with the same bin width (256 bins). 

In Figure Q] we show representative results obtained 
from MHD simulations with a range of Mach numbers and 
values of the plasma /3. We find that the PDF of 2:3 is recon- 
structed with rather good accuracy, even though the forms 
of the PDFs are only approximately lognormal, as we as- 
sumed in the calculation of £. The tails of the PDF (es- 
pecially the negative tails) are not well-recovered, but we 
note that A 3 data points are included in the calculation of 
the PDF of 111(0:3) while only A 2 data points are included 
in the calculation of the PDF of \n(x3R). Thus the extreme 
values of the field are less likely to be present in X3R than 
in X3, explaining the lack of sensitivity of the reconstructed 
PDF to the tails of the true PDF. The worst PDF recovery 
is for the M = 2, /3 = 0.1 simulation, which contains the 
highest degree of anisotropy (caused by the strong magnetic 
field, which produces sub-Alfvenic turbulence, as discussed 
in BFP). In general, the success of the PDF recovery is bet- 
ter at higher Mach number. 

Taking into account the understood lack of sensitivity 
to the PDF tails, the success of the PDF reconstruction im- 
plies that the PDF of 111(2:2) is, to a good approximation, 
a scaled, shifted copy of the PDF of 111(0:3) _ i.e. that the 
log-space PDF retains its form under projection, but is com- 
pressed due to line-of-sight averaging. BFP provide a means 
of determining the amount of compression (via R or, equiv- 
alently, £), which is dependent on the power spectrum of 2:3, 
and measureable (under the assumption of isotropy) using 
the power spectrum of x%. We find that, for the simula- 
tions analysed here, £ is typically around 2.7, but varies by 
about 0.5 around this value. Using an ensemble of slightly 
lower resolution numerical simulations (A = 128), driven at 
a range of spatial scales, Brunt & Mac Low (2004) found em- 
pirically that £ w 3. The A = 1024 simulations of Federrath 
et al (2009) yield £ w 2.9 and £ f» 2 for solenoidal and com- 
pressive forcing respectively. (In general, £ will depend on 
the scale ratio, A, and the form of the power spectrum of the 
density field, and should be calculated in each application.) 

The above results demonstrate that, to a good approx- 
imation, a lognormal density PDF projects into a lognor- 
mal column density PDF. However, one may rightly ques- 
tion the utility of the rescaling procedure. If we know the 
PDF is lognormal, and we have an estimate of the variance, 
it is straightforward to generate an analytic expression for 
the PDF. Indeed, we find that analytic lognormal functions 
scaled to a variance of ct^ are good representations of the 
density PDFs displayed in Figure [1] The true utility of the 
method must therefore be demonstrated for density PDFs 
that deviate significantly from a lognormal form. We con- 
duct this test in the next section. 

3.2 Self-Gravitating Hydrodynamic Simulation 

We now test the method on a density field produced by 
a simulation of star formation in self-gravitating hydrody- 
namic turbulence (Price & Bate 2009). Density field snap- 
shots were taken at 0.1, 0.2, 0.9, 1.0 free-fall times. At 
the earliest times during the development of the turbulence 
and also at t ~ tff once gravitational collapse has set in, the 
density PDFs differ substantially from a lognormal form and 
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Figure 1. Comparison of true 3D PDFs of ln(a;3) (dashed red lines) with reconstructed PDFs of ln(x3^j) (black lines) generated from 
information contained solely in the projected column density field for the MHD simulations. The parameters of each model (Mach 
number, M, and plasma /}) are quoted in each panel. 



therefore provide a good test of our method. The projected 
column density fields were analysed using the BFP method 
and estimates of the 3D normalised density variances were 
made (only one projection axis was used). We rescaled each 
X2 field by matching the variance of ax\ to cr^ 3i? , where £ is 
found by a series of sequentially refined test values. The re- 
sulting PDFs are shown in comparison to the true 3D PDFs 
in Figure [2] In this figure, we also show analytic lognormal 
functions with variance set to o~^ 3R for each field. For refer- 
ence, on each plot, we also draw vertical lines to mark the 
mean density and one-hundredth of the mean density. 

The evolution of the density PDF with time is as fol- 
lows. At early times (i < 0.6 tff) the PDF has an extended 
negative tail. A roughly lognormal form develops as t ap- 
proaches tff, accompanied by the development of an ex- 
tended positive tail. Such extended positive tails are seen in 
column density PDFs and are associated with gravitational 
collapse and star formation activity (Klessen 2000; Feder- 
rath et al 2008; Kainulainen et al 2009). 

We find that, in general, the density PDF is recovered 
very well at densities above the mean density. For most snap- 



shots, the rescaled column density PDF is a reasonably good 
representation of the density PDF at densities as low as one- 
hundredth of the mean density, but deviates from it substan- 
tially below this. The extreme low density tail is difficult to 
measure in most situations (especially in observational con- 
ditions which include noise) so this may not be a significant 
issue. For all snapshots (except t — 0.8 tff) the rescaling 
method performs better than using the analytic lognormal 
function. From the small number of snapshots available, it 
appears that the positive tail seen at late times is relatively 
more prominent in the rescaled column density field PDF 
than in the true density field PDF. This suggests that ob- 
served positive tails in column density PDFs imply the pres- 
ence of similar positive tails in density PDFs, but that some 
caution in their interpretation should be applied. 

The overall success of the method supports the idea that 
the form of the density PDF is preserved in the column den- 
sity PDF (even for a non-lognormal form) with appropriate 
provisos on the extreme positive and negative tails discussed 
above. 
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Figure 2. Comparison of reconstructed PDFs of ln(x3n) (black lines) with the true 3D PDFs of ln(z3) (red dashed lines) for the self- 
gravitating hydrodynamic simulation. The time measured in free-fall times is shown in each panel. The dotted blue lines show analytic 
lognormal functions with variance &x 3R f° r comparison. The vertical lines mark the mean density and one-hundredth of the mean density 
in each PDF. 



4 SUMMARY 

We have introduced and tested a simple method for recon- 
structing the probability density function (PDF) of a 3D 
turbulent density field using information present solely in 
the projected (observable) column density field in 2D. The 
method builds on a previously established method to calcu- 
late the 3D normalised density variance, recently presented 
by Brunt, Federrath, and Price (BFP, 2010). 

To a good approximation, the PDF of log(normalised 
column density) is a compressed, shifted version of the PDF 
of log(normalised density), but can deviate significantly in 
the extreme tails. The compression factor, £, can be derived 
observationally from the column density power spectrum, 
assuming statistical isotropy, using the BFP method. 
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